Instability of speckle patterns in random media 
with noninstantaneous Kerr nonlinearity 
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Onset of the instability of a multiple-scattering speckle pattern in a random medium with 
Kerr nonlinearity is significantly affected by the noninstantaneous character of the nonlinear 
medium response. The fundamental time scale of the spontaneous speckle dynamics beyond the 
instability threshold is set by the largest of times To and tnl, where Tb is the time required for 
the multiple-scattered waves to propagate through the random sample and tnl is the relaxation 
time of the nonlinearity. Inertial nature of the nonlinearity should complicate the experimental 
observation of the instability phenomenon. 
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Instabilities, self-oscillations and chaos are widely en- 
countered in nonlinear optical systemsi2*i In particu- 
lar, the so-called modulation instability (MI), consist- 
ing in the growth of small perturbations of the wave 
amplitude and phase, resulting from the combined ef- 
fect of nonlinearity and diffraction, can be observed for 
both coheren1i2i4 and incoheren1s£i§iL2ai light. In the lat- 
ter case, propagation of a powerful, spatially incoherent 
quasimonochromatic light beam through an optically ho- 
mogeneous nonlinear crystal results in ordered or disor- 
dered spatial patterns^iLSiS. Recently, instabilities have 
also been predictedASiil for diffuse light in random me- 
dia with Kerr nonlinearity (sample size L 3> mean free 
path £ ^> wavelength A). The strength of the nonlin- 
earity has to exceed a threshold for the instability to 
develop, similarly to the case of incoherent MI (coherent 
MI has no threshold), but in contrast to the latter, (a) 
initially coherent wave loses its spatial coherence inside 
the medium, due to the multiple scattering on hetero- 
geneities of the refractive index, (b) temporal coherence 
of the incident light is assumed to be perfect (in the case 
of incoherent MI, the coherence time of the incident beam 
is much shorter than the response time of the medium), 
(c) incident light beam is completely destroyed by the 
scattering after a distance ~ t and light propagation is 
diffusive in the bulk of the sample, and (d) scattering 
provides a distributed feedback mechanism, absent in the 
case of MI. 

Since the noninstantaneous nature of nonlinearity 
proved to be important for the development of optical 
instabilitiesiSiSiiiSiSiLSiS, it is natural to consider its ef- 
fect on the instability of diffuse waves in random non- 
linear media. Up to now, the latter phenomenon has 
only been studied assuming the instantaneous medium 
respons e) 1 ?' 11 . In the present letter, we show that while 
the absolute instability threshold does not depend on 
the nonlinearity response time tnl, the dynamics of the 
speckle pattern above the threshold is extremely sensi- 
tive to the relation between t nl and the time To = L 2 / D 
that takes the diffuse wave to propagate through the ran- 



dom sample {D — ci/3 is the diffusion constant and c is 
the speed of light in the average medium). We derive 
explicit expressions for the characteristic time scale r of 
spontaneous intensity fluctuations beyond the instabil- 
ity threshold and the maximal Lyapunov exponent A max , 
characterizing the predictability of the system behavior. 
Finally, we discuss the experimental implications of our 
results. 

We consider a plane monochromatic wave of frequency 
too incident on a random sample of typical size L £. 
The real dielectric constant of the sample can be written 
as e(r, t) = eo + Se(r) + A£nl(i", t), where eo and 5e(r) are 
the average value and the fluctuating part of the linear di- 
electric constant, respectively, and A£nl(i", t) denotes the 
nonlinear part of e. Without loss of generality, we assume 
eo = 1 and (&(r)&(ri)} = 47r/(fco^)<5(r— ri). Under con- 
ditions of weak scattering kg£ 3> 1 (where fco = uq/c), 
the average intensity (I(r)) of the scattered wave inside 
the sample can be described by a diffusion equation, the 
short-range correlation function of intensity fluctuations 
dl(r,t) = I(r,t) — (I(r)) decreases to zero at distances 
of order A, and the long-range correlation of SI(r, t) can 
be found from the Langevin equation^ (see Ref. [Kj 
a review of wave diffusion in random media) : 



for 



St 



5I(t, t) - DV 2 SI(r, t) = -V ■ j ext (r, t). (1) 



In a nonlinear medium, the random external Langevin 
currents j ex t (r, t) obey a dynamic equationAi 

d [ f°° 

^JcxtM) = y y d3r 'y n dAiq(r,r',At) 



x — Ae NL (r',t- At), 



(2) 



where the integration is over the volume V of the random 
medium and q(r, r', At) is a random, sample-specific re- 
sponse function with zero average and a correlation func- 
tion given by the diagrams of Fig. ^ 

Supplemented by a Debye-relaxation relation for 
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FIG. 1: Diagrams contributing to the correlation function 



(r, r ', t — t ')q ,(j) * (ri , r{ , t\ — t[)j of the random response 

functions q entering into Eq. 0. Solid lines represent the 
retarded and advanced Green's functions of the linear wave 
equation, dashed lines denote scattering of the two connected 
wave fields on the same heterogeneity. Wavy lines are k$ 
vertices. 
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FIG. 2: Instability threshold as a function of the excitation 
frequency O (in units of TpJ 1 ) for several ratios of the non- 
linearity relaxation time tnl and the time Tb required for 
light to diffuse through the random sample. Dashed line cor- 
responds to tnl /To = 0. Inset: Lyapunov exponent A vs. Q 
(both in units of T£ ) slightly above the threshold (p = 1.1) 
for the same values of tnl /To as in the main plot. 



Ae N L(r,£): 
d 

r NL — Ae NL (r, t) = -Ae NL (r, t) + 2n 2 J(r, i), (3) 
at 

where ri2 is the nonlinear coefficient, Eqs. and 
can now be used to perform a linear stability analysis 
of the multiple-scattering speckle pattern in a random 
medium with noninstantaneous nonlinearity. Assuming 
SI(r, t) = SI(r, a) exp(crf) with a = iCl + A (Q > 0), and 
similarly for j ex t( r j£) and Ae NL (r, t), we use Eqs. 
and J2Jl to express the correlation function of Langevin 

currents /iext( r > a )iext* ( r i > a f) m terms of the intensity 
correlation function (6I(r,a)6I*(ri,a)) and Eq. to 
express the latter correlation function in terms of the 
former. For large enough sample size L/£ ^> fco^ and 
moderate frequency flT® <C [L/(fco^ 2 )] 2 we can neglect 
the short-range correlation of intensity fluctuations and 
write the condition of consistency of the two obtained 
equations [and hence the condition of consistency of Eqs. 

mm] as 

P ~F(nT D ,AT D )iT(fiT NL ,AT NI ^ (4) 

where a numerical factor of order unity has been omit- 
ted, p = (Att-nl) 2 (L/£) 3 is the effective nonlinearity pa- 
rameter, (Atinl) = "-2 (I( r )) is the average value of the 
nonlinear part of refractive index, and (/(r)) is assumed 
to be approximately constant inside the sample. The 
function F in Eq. (@J is the same as in the case of in- 
stantaneous nonlinearity 11 , while H (x, y) = x 2 + (1 + y) 2 
originates from the noninstantaneous nature of the non- 
linear response. In the limit of t^l, (instantaneous 
nonlinearity), H = 1 and Eq. Q reduces to Eq. (7) of 
Ref.m 



For given f2 and p, Eq. Q determines the correspond- 
ing value of the Lyapunov exponent A. A > signifies 
the instability of the speckle pattern with respect to the 
excitations at frequency il. If p < 1, then A < for 
all f2. As p exceeds 1, A becomes positive inside some 
frequency band (0,r2 max ) (see the inset of Fig. |2J). We 
show the frequency-dependent threshold following from 
Eq. igj in Fig. for different values of t nl /T d . The 
threshold value of p at > t^"l deviates from the result 
obtained for the instantaneous nonlinearity (dashed line 
in Fig. |2J), while the absolute threshold is p = 1 at £1 = 0, 
independent of tnl/^d- 

To study the onset of the instability, we consider the 
limit of < p — 1 <C 1 in more detail. In this limit, 
a series expansion of the function F in Eq. Q can be 
found: F(x,y) ~ 1 + C\x 2 + y(C 2 — C 3 x 2 ), where we 
keep only the terms quadratic in x and linear in y, and 
the numerical constants are C\ ~ 0.031, C2 — 0.59, 
C3 ~ 0.045. At a given frequency O, the threshold value 
of p then becomes p ~ 1 + [d + (t nl /T d ) 2 ] (^T d ) 2 , 
while at a given p the maximum excited frequency is 
= T^ 1 [Ci_+(t nl /T d ) 2 ]" 1/2 (p_- 1)V2 and the 
shortest typical time of spontaneous intensity fluctua- 
tions is t — l/f2 max . 

As follows from the above analysis, a continuous low- 
frequency spectrum of frequencies (0,f2 max ) is excited at 
p > 1, and the Lyapunov exponent A decreases mono- 
tonically with (see the inset of Fig. This allows 
us to hypothesize that at p = 1 the speckle pattern un- 
dergoes a transition from a stationary to chaotic state. 
Such a behavior should be contrasted from the "route to 
chaos" through a sequence of bifurcations, characteris- 
tic of many nonlinear system a 1 ' 2 ' 3 . Given initial condi- 
tions at t — and deterministic dynamic equations, the 
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evolution of a chaotic system can only be predicted for 
t < 1/A max . We find the maximal Lyapunov exponent 
A max = Tq 1 (p-1)/(C 2 + 2t nl /T :d ). It is now easy to see 
that the fundamental time scale of the speckle pattern dy- 
namics at p > 1 is set by the largest of the times Td and 
tnl- Indeed, in the case of fast nonlincarity (tnl "C lb) 
we have ft max oc T I) 1 (p - l) 1 / 2 , r oc T D (p - 1)~ 1/2 , and 
A m ax oc Tp 1 ^ — 1). In the opposite limit of slow nonlin- 
ear response (tnl 3> 2b) we find f2 max oc Tj7l(p — l) 1,/2 5 
t oc tnl(p- 1)~ 1/2 , and A max oc T^(p — 1), respectively. 

The origin of the instability of the speckle pattern 
J(r, t) in a nonlinear random medium and its sensitivity 
to tnl can be qualitatively understood by considering / 
as a result of interference of many partial waves traveling 
along various diffusion paths. An infinitely small pertur- 
bation of I at some point ro at time to diffuses throughout 
the medium [according to Eq. l|T|l] and, after some time 
At, modifies the speckle pattern inside a sphere of radius 
R ~ (DAt) 1 / 2 around ro. Due to the nonlinearity of the 
medium, this leads to a change of the nonlinear part of 
the dielectric constant Aenl inside the same sphere after 
a time ~ tnl [according to Eq. ©] and, consequently, to 
a modification of the phases of the partial diffuse waves 
that interfere to produce a new speckle pattern I(r,t). 
If the strength of the nonlinearity exceeds the threshold 
given by Eq. Q, I(r ,t) — I(r 0} t ) can be larger than 
the initial perturbation. The latter is hence amplified 
and the speckle pattern develops spontaneous dynam- 
ics. Noninstantaneous nature of the nonlinearity leads 
to the damping of the high-frequency components in the 
spectrum of the initial intensity perturbation, raising the 



instability threshold for > t^t , as we show in Fig. [2] 
This phenomenon is similar in its origin to the increase of 
the speed of coherent MI patterns with decreasing tnl m 
a homogeneous nonlinear medium*. Similarly to the case 
of coherent Mis, the instability is arrested in the limit 
of tnl — * oo, in contrast to the case of incoherent MI, 
where tnl — > oo (more precisely, tnl 3> than the coher- 
ence time of the incident beam) is an explicit assumption 
of the theoretical model o 5 i 6 i 9 . 

To conclude, we discuss the experimental implications 
of our results. Atinl up to 10 -3 can be realized in in- 
organic photorefractive crystalsS*^. If scattering centers 
(optical defects) are introduced in the crystal, Ljl ~ 10 2 
will suffice to reach the instability threshold p ~ 1. Al- 
ternatively, in nematic liquid crystals both nonlinear- 
ity (Atinl up to 0.1, see, e.g., Ref. H3) and scattering 
(L/£ > 10 has been reported in Ref. lint ) can be suffi- 
ciently strong to obtain p ~ 1. The nonlinearity is rather 
slow in the above cases, tnl exceeds Td (7b ~ 10 -j- 100 
ns for t ~ 0.1 -j- 1 mm and L ~ 1 cm), and the results 
obtained in the present letter are, therefore, of partic- 
ular relevance. Finally, although the noninstantaneous 
nature of nonlinearity does not affect the absolute insta- 
bility threshold, it should complicate the experimental 
observation of the instability phenomenon. Indeed, in a 
real experiment only sufficiently fast fluctuations can be 
detected. Inasmuch as i7 max decreases with tnl, larger 
p (and thus stronger nonlinearities) will be required for 
the instability to be detectable in a medium with nonin- 
stantaneous nonlinearity. 
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